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We discuss quantum entanglement in the context of the thermodynamic arrow of time. We 
review the role of correlations in entropy-decreasing events and prove that the occurrence of a 
transformation between two thermodynamic states constitutes a new type of entanglement witness, 
one not defined as a separating plane in state space between separable and entangled states, but as a 
physical process dependent on the local initial properties of the states. Extending work by Partovi, 
we consider a general entangled multipartite system that allows large reversals of the thermodynamic 
arrow of time. We describe a hierarchy of arrows that arises from the different correlations allowed 
in a quantum state and examine these features in the context of Maxwell's Demon. We examine in 
detail the case of three qubits, and also propose some simple experimental demonstrations possible 
with small numbers of qubits. 
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INTRODUCTION 



For over a century, Maxwell's demon has provided a 
setting in which to address the limitations that thermo- 
dynamics places on an observer free to perform mea- 
surements on a system and then act on their acquired 
information in some algorithmic way As is now 

well appreciated, information has an energetic value. A 
demon Q in possession of information about a physi- 
cal state may transform this information into mechanical 
work [H| . Specifically, given a system of dimension D 
in a state p and the presence of a reservoir at a temper- 
ature T, the demon can extract W = kT (hi D — S[p\) 
amount of mechanical work from the reservoir, where 
S[p) = — Tr[plogp] is the von Neumann entropy of the 
state. The demon fails in his attempt to violate the sec- 
ond law of thermodynamics because the demon's mem- 
ory, where he records the measurement results of the 
state, must be erased in order to operate in a cycle [6j. 
Such an erasure of memory can only occur with an un- 
avoidable dissipation of heat Q . 

As a concrete example, we may consider a gas of N 
particles that resides in a piston on one side of a mov- 
able partition and in thermal equilibrium with a reservoir 
at a temperature T, Fig. [TJ The demon, upon discover- 
ing what side of the partition the gas is on, may extract 
usable energy by deterministically acting on the piston 
system - if the gas is to the left of the partition, the de- 
mon puts the piston in contact with the reservoir, loads 
the partition with a mass M, tilts the piston as in Fig. [T] 
and allows the gas to expand slowly through the available 
empty volume. The expansion of the gas does work on 
the mass, elevating to a height h, which makes mechan- 
ical energy available. Indeed the expansion can provide 
a maximum mechanical energy of 7V/cTln2, which is ex- 
tracted from the thermal reservoir. 



Confounding MaxwelPs Demon 

In order to illustrate the new possibilities that quan- 
tum entanglement can bring to the story, it is amusing 
to consider cases in which a "global demon" , with access 
to a large entangled quantum state, can confuse the tra- 
ditional "local demon" who can only measure and act lo- 
cally. The global demon can arrange that any thermody- 
namic process for the local demon can run "backwards" 
- heat can flow from a cold to hot and gases can contract 
instead of expand - driven by entanglement present in 
the state. 

As a concrete example, consider a system that consists 
of a reservoir and piston that are both in local thermal 
states, but with the global state of the system being pure. 
This state is prepared by the global demon, who may 
manipulate the global state as he pleases. The global 
demon ensures that the entanglement is local, in the sense 
that any part of the total system is entangled with its 
surroundings. Furthermore, for simplicity, we constrain 
the reservoir subsystems to couple to each other much 
weaker than to the piston subsystem [35| . This situation 
is indicated in Fig. [TJ 

Since thermality arises from entanglement between the 
different components of the system in a global pure state, 
when the local demon places the piston in contact with 
the quantum reservoir and releases the partition, instead 
of the gas expanding up the piston to fill the empty vol- 
ume, the thermal state repeatedly expands forwards and 
then contracts backwards in the tube as energy flows into 
and out of the piston. The spontaneous contraction phase 
of the gas is identical to the usual spontaneous expansion 
of a gas but with time "running backwards" . 

For the local demon, the thermodynamic arrow of time 
has failed and he is unable to extract any net work despite 
having been given a piston that is on its own completely 
indistinguishable from a standard thermal state on one 
side of a partition, seemingly coupled to a standard ther- 
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FIG. 1: Uncorrelated case: The piston is in contact with a 
thermal reservoir of particles, at temperature T. Mechanical 
energy is extracted from the reservoir by the demon choosing 
which way to tilt the piston, conditional on his information 
about what side the gas is located. 




FIG. 2: Entangled case: The local thermodynamic states 
of the piston and reservoir are indistinguishable from the un- 
correlated cases, however now entanglement is present in the 
system. In this case the piston both expands forwards and 
contracts backwards. Heat flows both into and out of the 
piston and the usual thermodynamic arrow no longer holds. 



problem of resolving the thermodynamic arrow's origin 
requires us to account for the particular initial conditions 
of our universe [8j . 

The notion of an "arrow of time" is quite broad - ex- 
tending from the uniform large-scale expansion of the 
universe, to the radiative arrow, and through to the indi- 
vidual, psychological perception of time [9|,[lQ|. However, 
since the laws of physics are invariant under time rever- 
sal (strictly CPT invariant), it is believed that any form 
of temporal directionality will stem ultimately from the 
initial conditions of the universe. 

Recent work by Partovi [TBI has considered a specific 
pure entangled state that produced heat flow from a cold 
body to a hotter one, contrary to the usual thermody- 
namic arrow. In [TBI Partovi identifies the large degree 
of entanglement present in the state as being responsi- 
ble for this unusual behaviour. We will see, in fact, that 
the negative heat flow he considers does not necessar- 
ily depend on the presence of entanglement in the state 
and, by adopting an information theoretic framework, we 
explore these and similar possibilities as a way of eluci- 
dating the properties of entangled quantum states and 
the connection with varied notions of "arrows of time" . 



mal reservoir. 

The mischievous global demon, however, may further 
confuse the local demon by thermalizing the entangle- 
ment in the (quantum) reservoir that he has access to, 
using a different (fully thermal, classical) reservoir. Since 
the full system of dimension D is in a pure state, the 
global demon can extract kThiD amount of work from 
the classical reservoir, which he may store as mechani- 
cal energy, leaving the quantum reservoir plus piston in 
a maximally mixed state. Then using only some of this 
mechanical energy he may convert the maximally mixed 
reservoir into a true thermal state at the original temper- 
ature T. At which point the thermodynamic arrow is re- 
turned, the global demon has converted the entanglement 
into mechanical energy and the local demon, unaware of 
all the entanglement that was present, will finally be able 
to extract NkT In 2 of work from the reservoir as he had 
originally hoped. 

The Arrow of Time 

There are many reasons to be interested in the arrow 
of time beyond simply fun stories about demons. The 
standard thermodynamic arrow applies to systems that 
are in isolation and amounts to a directionality for cer- 
tain physical processes. The curious fact that certain 
processes are forbidden and certain states inaccessible is 
generally traced back to the initial conditions of the phys- 
ical system. Furthermore, since all parts of the observed 
universe seem to obey a common "arrow of time", the 



The Thermodynamic Arrow and Correlations 

A paradigmatic setting that captures the essentials of 
the thermodynamic arrow of time is in the context of 
heat flow between two systems A and B that interact 
in isolation, so that their total energy is constant. For 
simplicity, we assume that they do no mechanical work on 
each other, but may exchange energy in the form of heat, 
for which conservation of energy implies that Qa + Qb = 
0, where Qa is the heat acquired by A and Qb is the heat 
acquired by B. 

Ultimately, the key ingredient that dictates which 
states are accessible via the set of all interactions on 
the total state pab is the degree of correlations ini- 
tially present between A and B. These correlations 
between the two subsystems are quantified by the mu- 
tual information, which is non-negative for all states 
and defined in terms of relative entropy as I(A : B) = 
S[pab\\pa ® Pb] = S A + S B - Sab, where, for exam- 
ple, Sa = — Tr[p^logp^] is the von Neumann entropy 
for the state pA- The mutual information is a convenient 
measure of how distinguishable the state pAB is from the 
completely uncorrelated product state pa® Pb- 

Of particular interest is when the reduced states at 
the initial time t{ of subsystems A and B are thermal 
states for their respective individual Hamiltonians, Ha 
and Hb- Given temperatures Ta and Te, the reduced 
states then take the form pa(U) oc exp[— (3aHa\ and 
Pb{U) oc exp[— PbHb], where fix = 1/kTx for subsys- 
tem X e {A, B}. 

We may consider the most general process of switching 
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on time-dependent interactions Vab (t) between A and B 
that evolve the composite system into a new state pABif)- 
For subsystem X G {A, B} at a later time f, the two 
thermodynamic variables of central interest are the av- 
erage energy Ux — ^[Hxpxif)} and the entropy Sx — 
—Tr[px(t) log px(t)]. Locally at A the interaction results 
in the transformation pa(U) pA(tf) = TYb[pab (£/)]> 
and so U a -» C^a + and SU ->> SU + ASU, with 

similar expressions for B. 

These changes in local thermodynamic variables 
are constrained by the fact that a thermal state 
exp[— PH]/Tr[ex.p[— f3H]] minimizes the free energy func- 
tion F[p] = Tr[pH] — S\p]//3. Consequently we have that 
Tx[p(U)H] - S[p(U)]/f3 < Tx[p(t f )H] - S\p(t f )]/P, which 
more explicitly yields 



PaAUa-ASa > 
(3bAUb-ASb > 



(1) 



for the change in the thermodynamic variables defined at 
A and B. It must be emphasized that equations (pQ) are 
kinematical restrictions on any transformation from ini- 
tially thermal states, {pa(U), Pb(U)}, to any other pair 
of local states {/?a (£/"), PB(tf)}- 

We consider transformations on the composite system 
that satisfy S[pab(U)] = S[pAB{tf)) (which includes uni- 
tary transformations) and so AI(A : B) = AS a + ASb- 
The only other restriction that we place on the transfor- 
mation is that Ty[pab(U)(Ha+H b )} = Ti[pAB(tf)(HA + 
Hb)]- Since we don't consider mechanical work, this 
amounts to requiring that there is zero net heat flow into 
the composite system, and the total energy before the 
interaction process equals the total energy afterwards. 
Writing Qa = AUa and Qb — AUb for the heat gained 
by A and B respectively, we obtain for these energy pre- 
serving processes that 



PaQa + PbQb > A/(A : B) 



(2) 



where AI(A : B) is the change in the mutual information 
between A and B. 

The special initial state pab{U) — Pa(U) ® Pb(U) i n 
which correlations vanish was historically called the con- 
dition of Stosszahlansatz, or 'molecular chaos'. If we con- 
sider A and B initially uncorrelated at some time U then 
I(A : B\ti) = 0, however interactions allow correlations 
to form between the two subsystems, and so at later times 
tf we have I (A : B\ tf) > 0. We immediately deduce that 
PaQa + PbQb > 0, or more explicitly in terms of tem- 
perature 
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Tb 



> 0. 



(3) 



In other words, heat can only flow from hot thermal 
states to cold thermal states in isolation. This is the 
standard thermodynamic arrow of time. 



This situation can be generalized to N initially uncor- 
related thermal states and provides us with the constraint 
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Qj 



> o 



(4) 



which forbids certain types of heat flow; for example, if 
we partition the TV thermal states in two, such that all 
the temperatures in the first group are lower than those 
in the second, then it is impossible for heat to flow from 
the first to the second. 

In general, however, initial correlations are expected 
to be present, I(A : B\U) > 0, and so the change in 
mutual information can be negative. From (|2]) it is clear 
that there is no longer a constraint on the directional- 
ity of heat flow between A and B [16]. As we shall now 
see, the correlations that make up the total mutual in- 
formation can arise from classical correlations, or from 
a combination of classical correlations and quantum en- 
tanglement. 



Entanglement correlations 

Entangled quantum systems can possess far stronger 
correlations than are possible classically, and in fac t the 
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vast majority of quantum states are entangled 
Work has been done previously on how quantum cor- 
relations affect macroscopic thermodynamic observables 
such as the susceptibility or the heat capacity, and it 
has been shown that measurement of these macroscopic 
observables can act as entanglement witnesses [13, fig . 
We shall show that the thermodynamic transformations 
themselves also constitute entanglement witnesses. We 
should be careful to distinguish this research p rog ram 
from the one to do with with formal analogies [12|, LL4( be- 
tween irreversible entanglement transformations and the 
second law of thermodynamics, on which much progress 
has recently been made [ll| , or from attempts to connect 
the second law of thermodynamics to Bell inequalities 
[l3( | . Here we are interested in the physical effects which 
entanglement between systems and reservoirs has on the 
thermodynamical transformations of the systems. 

Recently in [l5| Partovi considered how the presence of 
entanglement in a system affects irreversible thermody- 
namic transformations. The intriguing scenario he con- 
sidered involves two subsystems A and B that are over- 
all in a pure state \^ab) and possess strong quantum 
correlations. Furthermore, he arranged that the local 
reduced states Pa,b — ^b,a[\^ ab}(^ ab\] are thermal 
states for the two subsystems and so locally indistin- 
guishable from classical thermal mixtures. Despite be- 
ing locally indistinguishable from classical thermal mix- 
tures, the composite system behaves quite differently un- 
der any energy-conserving unitary. Since the total ini- 
tial state is a pure state we have that S[pab{U)] — 0, 
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which remains true at any later time. From a Schmidt 
decomposition for A and B we find that pA(t) and pB(t) 
are isospectral, Spec(pA^)) = Spec(/0e(£)), and so their 
entropies are always equal. Consequently, any unitary 
transformation occurring on the composite system will 
have AS a = ASb- 

Once again, assuming (i) no mechanical work and (ii) 
no overall change in the total energy, the condition on the 
heat flow between A and B is PaQa + PbQb > AS a + 
ASb- However, since AS a = ASb and Qa = —Qb we 
have that 



arrow in the form of energy transfer from a colder sub- 
system to a hotter one. 

The total correlations, as quantified by I (A : £>), also 
contain classical correlations, which may be defined as 



■AS a (1/Pa + 1/Pb) >0 



(5) 



and so AS a = ASb < 0. Energy conserving unitary 
interactions never increase the local entropies of either 
subsystem, in stark contrast to the uncorrelated situation 
of classically mixed thermal states, where ASa + ASb 
never decreases. 

However, the model in [15| does not actually require 
quantum entanglement. For example we may simply de- 
phase \^ab) to turn it into a classically correlated sep- 
arable state that would produce the exact same reversal 
of the thermodynamic arrow. In addition, while the pure 
state I^ab) allows local entropy decreases, and the re- 
versal of the thermodynamic arrow, the isospectral con- 
straint for the bipartite splitting of the pure state \^ab) 
is very restrictive - the only way in which two different 
temperatures for A and B is achieved in [15| is to impose 
that the energy spectrum of A be a scaled version of the 
energy spectrum of B. Another undesirable feature of the 
model is that the initial thermal states that are consid- 
ered in [l5| (squeezed Gaussian states of two oscillators) 
are the only thermal states that are attainable in this 
setting - all other thermal states are inaccessible. 

In the next section we extend the work of [TBI to show 
that we may distinguish entangled systems from classi- 
cally correlated system through violations of thermody- 
namic arrow and produce transformations between arbi- 
trary thermal states without restrictions on the energy 
spectra. Furthermore, we analyze quantum systems in 
which the correlation structure is more complex and al- 
lows for a more subtle range of violations of the thermo- 
dynamic arrow. 



Classical correlations and entanglement witnesses 

Recently [lo[ the mutual information I(A : B) was 
given an operational meaning in terms of the minimal 
amount of local randomizing work that must be done 
on a bipartite system to reduce it to a product state 
and destroy all correlations. We will give here another 
thermodynamical interpretation - namely, the mutual in- 
formation is a measure of the maximal amount that the 
isolated bipartite system may violate the thermodynamic 



I C (A : B) = m&K MA ®M B H(A : B) 



(6) 



where H(A : B) is the classical (Shannon) mutual infor- 
mation for the joint probability distributions generated 
by local POVM measurements M A ®M B at A and B (22| . 
The classical mutual information I C (A : B) is always less 
than or equal to the total mutual information I(A : B) 
and vanishes if and only if the bipartite state is a product 
state [lal . 

It is clear that for bipartite states pab the classical 
mutual information always obeys the bound < I C {A : 
B) < logD, where D is the dimensionality of the smaller 
subsystem. On the other hand, while the quantum mu- 
tual information I(A : B) is always positive, it can in 
general take on values larger than log!}. However, if we 
write the quantum mutual information as 

I(A :B) = Sa + Sb- Sab 

= Sa — (Sab — Sb) = Sb — (Sab — Sa) 

and restrict ourselves to separable states pAB for which 
Sab is always greater than both Sb and Sa [23j we 
may deduce that I (A : B) < min{SU, Sb} and so 
I (A : B)\ sep < \ogD over the set of separable states. This 
upper bound is saturated for perfectly correlated, zero 
discord, separable states pAB — tj J2i K)a(^| ® in 
which case I C (A : B) = I (A : B) = \ogD. 

However, for entangled states we can have I (A : B) > 
logD, and so from (|2]) we see that certain types of heat 
flow are possible that are forbidden classically. In par- 
ticular a quantity of heat I (A : B)/(\/3a — Pb\) can be 
transferred from the cold subsystem to the hot system, 
whereas classical correlations could only permit at most 
a quantity of heat \og D / (\f3 a — Pb\) to flow from cold 
to hot. From this it follows that the reversal of the ther- 
modynamic arrow of time is a new kind of entanglement 
witness where quantum correlations present in the state 
allow much larger violations of the arrow than are classi- 
cally possible, and so the detection of such large reversals 
must imply an entangled state. Standard entanglement 
witnesses correspond to an observable W that defines a 
plane Tr[Wp] = in state space separating entangled 
states from separable ones, while for us, it is the set of 
energetically accessible states that acts as the witness. 
Instead of local observables, it is the transformations of 
the state that signal entanglement. 

The classical mutual information I C (A : B) has pre- 
viously been considered in the context of the thermo- 
dynamic arrow as a candidate measure for the classical 
memory record that a system A has of events that affect 
B [27j| . The claim in [27] was that any event that de- 
creases the entropy of B necessarily reduces the classical 
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mutual information between A and and so entropy- 
decreasing events at B do not leave a memory trace in 
A. This was a proposed resolution to the empirical fact 
that we only observe entropy-decreasing events, namely 
the "arrow of time dilemma" . 
However, 



m 



28| we demonstrate explicitly that this 



claim is false, and in reality quantum mechanics allows 
the exact opposite to occur. We show that large quan- 
tum correlations present in a state can be used up to re- 
verse the thermodynamic arrow and at the same time in- 
crease the classical correlations between two subsystems. 
Quantum mechanics allows enhanced memory records of 
entropy-increasing events and so there can be no resolu- 
tion of the arrow of time in terms of classical memory 
records. 

The process of dividing correlations into quantum 
and classical components has many subtle aspects and 
counter-intuitive results exist, for example in the setting 
of multipartite correlations it is possible to have purely 
quantum multipartite correlations without any classical 
multipartite correlations being present [26|. In the next 
section we consider multipartite quantum systems that 
are in a highly entangled pure state and analyse how the 
quantum correlations present allow arbitrarily large re- 
versals of the thermodynamic arrow. 



LOCAL THERMAL STATES, GLOBAL STATES 
AND A HIERARCHY OF ARROWS 



We now turn to an analysis of entanglement-assisted 
violations of the thermodynamic arrow in a general mul- 
tipartite setting. Once again, we consider an isolated 
system that may undergo a complex time-dependent in- 
teraction process between times U and tf, but with the 
total energy, defined via a total Hamiltonian iJ tot , at time 
U equal to the total energy at time tf. 

It is well known that most quantum states are highly 
entangled 25], however recently a deep theorem Q has 
shown that if the total system is large, then any randomly 
choosen pure state |$) that satisfies the total energy con- 
straint has the property that any sufficiently small sub- 
system is highly likely to be in a thermal state. For 
generic large systems in a pure state, thermality of its 
subsystems naturally arises from the high degree of en- 
tanglement present in the state. Consequently, instead 
of entangled thermal states such as in [l5| being artificial 
or exceptional, for large systems in a randomly chosen 
pure state they are actually quite typical. 



Local thermal states from an entangled state 

For simplicity we shall consider the multiple subsys- 
tems Si, . . . Sjv, to be qubits, but it is clear that similar 
arguments apply to any higher dimensional subsystems. 



We take the reduced states, obtained by tracing out the 
other subsystems, to be initially thermal and the total 
system is assumed to be in a pure state \^f) with a fixed 
energy. We refer to the reduced states pi, P2, • • • as the 
thermal marginals for the subsystems. 

However, not all sets of thermal marginals may be ob- 
tained from a global pure state |^). For the case of 
qubits, if the marginal states for the individual qubits 
are given by {pi}, then the necessary and sufficient con- 
ditions for the existence of a pure state such that 
Tij^i\^/}(^/\ = pi are in terms of the smallest eigenvalue 
of each subsystem, = min(Spec(pi)). The conditions 
are [30|l 



< A; < 1/2 i 



1,...,N. 



(7) 



Without loss of generality, we take the Hamiltonians 
of each qubit to be equal, and choose the ground state 
to have zero energy while the excited state has energy 
1. Thus, Hi = |(J + Zi), where Zi is the Pauli op- 
erator for the Z direction, and for qubits with thermal 
marginals, the total expected energy of the system is then 
E = ^ZiXi- The parameter space for the set of states 
in terms of the smallest eigenvalues of its subsystems 
is defined by (0) and is an TV-dimensional polyhedron 
Vn-, with the constant energy condition being a hyper- 
plane that intersects Vn on a subset T/v-i of dimension 
at most N — 1. Each point in 7/v-i corresponds to an ac- 
cessible combination of thermal marginals, however there 
is in general more than one pure global state associated 
to such a point. 

As an example, in the case N = 3, the region of pa- 
rameter space is a diamond formed from two tetrahe- 
dra, while the constant energy condition corresponds to a 
plane that slices the diamond in triangular cross-sections, 
as depicted in Fig. 03 Each point in the diamond de- 
fines three qubit thermal states, however, since the over- 
all state is a pure state, any other pure state may be 
reached by a general unitary transformation. It is the re- 
striction to energy-conserving unitaries which means that 
only thermal states with parameters in the triangular re- 
gion ?2 are accessible. Such energy conserving unitaries 
could arise from time-dependent interactions Vij(t), or, 
for example, from time-independent interactions of the 
form XiYj -Y { Xj. 

The diamond polyhedron has some nice features. The 
origin is the bottom vertex and corresponds to E = 0, 
where all the qubits are at T = 0, while E = 3/2 is 
the other extreme point, where each qubit is maximally 
mixed, and T = oo. The centroid of any given constant 
energy triangle corresponds to all three qubits having 
equal temperatures and so would correspond to the stan- 
dard equilibrium one would expect if the three subsys- 
tems were initially uncorrelated. For standard thermally 
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FIG. 3: The parameter space V?,: The set of entangled 
states giving local thermal states possesses a diamond poly- 
gon parameter space. A triangular slice T2 with constant total 
energy is shown. Dark regions have low energies and temper- 
atures, while bright regions have high energies and tempera- 
tures. 
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FIG. 4: The parameter space 7*2: The standard thermo- 
dynamic arrow would constrain C to evolve only into the 
cone of states depicted, and this drives it steadily towards 
equilibrium, while the entangled state Q has no such re- 
striction. Dotted lines are the isotherms for T a , Tb and T c , 
while the solid internal lines represent (T a = Tb,T c = 0), 
(T b = T c ,T a = 0) and (T c = T a ,T b = 0). The center point 
corresponds to thermal equilibrium T a — Tb — T c . 



mixed states free to interact, the total system would tend 
to evolve towards this configuration. However, from our 
analysis we see that this is no longer imposed. The over- 
all state is pure and the total system can undergo unitary 
evolution to any other set of local thermal states in 75- 
The large degree of entanglement present has radically 
lifted the constraints of the thermodynamic arrow and 
allows the subsystems to transformation to otherwise in- 
accessible states. 

The triangle formed by the intersection of E = 1 and 
the diamond forms the widest part of the polyhedron and 
is special in the sense that there is enough entanglement 
for two of the qubits to form a singlet state, but not 
too much that the singlet becomes impossible. In terms 
of temperature, for points on this triangle it is possible 
to unitarily transform along E =constant to a situation 
where two qubits are maximally mixed while the third is 
in a pure state. For E < 1 it is always possible to turn 
off the temperature of one qubit with the other two being 
at some finite temperatures, while for E > 1 this is no 
longer possible and each qubit always has some non-zero 
temperature. 

The diamond setting also provides a simple way to 
visualise the contrast between a system subject to the 
standard thermodynamic arrow and the highly entangled 
system discussed here. Fig. |4] shows the constant en- 
ergy region T2 for a system of three qubits with thermal 
marginals p a , p^ and p c . The centroid of the triangle cor- 
responds to the equilibrium configuration T a = X5 = T c . 
The solid lines passing through the equilibrium point and 
intersecting the midpoint of the sides correspond to the 
situation where two of the subsystems have equal tem- 
peratures. These lines divide the triangle into six regions, 



labelled I- VI in Fig. [5j corresponding to the six choices of 
temperature orderings, T a < T5 < T c , • • • , T c < T5 < T a . 

If the system is in a thermal configuration C with no 
correlations then equation (|4]) forces the system to evolve 
only in the direction shown in Fig. HJ The thermody- 
namic arrow therefore imposes a cone of accessible con- 
figurations for the thermal system, and gradually drives 
the system inexorably to thermal equilibrium at the cen- 
ter of the region (a fun comparison is with black holes 
- the light cone of a person inside the event horizon al- 
ways forces them to move radially inwards, eventually 
reaching the central singularity, while here we have in 
some sense a 'thermal cone' that dictates how the sys- 
tem must evolve in time, ultimately reaching the central 
equilibrium point). 

It makes physical sense in this context to only con- 
sider a small transfer of energy that perturbs the initial 
thermal configuration into a neighbouring thermal con- 
figuration. A state with T5 = T c has no formal restriction 
from (|4j) on the amount of heat that may be exchanged 
between b and c, however the configuration T5 = T c is 
stable in the sense that a small transfer of heat will auto- 
matically induce a correcting directionality that returns 
the system to X5 =T C . 

This is illustrated in Fig. [5] where given an initial ther- 
mal configuration C, subject to random heat exchanges 
and the constraint (HJ) , the system starts in region I where 
T a < T5 < T c . Heat is exchanged, cooling c and heating 
a and b until the system crosses the X5 = T c line into 
region II. At which point only fluctuations around this 
line are permitted by (|4|). Subsequently, the system a 
is heated up until the thermal configuration reaches the 
central equilibrium point, as shown in the diagram, at 
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FIG. 5: Contrast in system evolution due to random trans- 
fers of energy between subsystems. For C, subject to the con- 
straint @ , the thermodynamic arrow forces the subsystems to 
the equilibrium configuration. The system Q is not restricted 
for the small fluctuations and performs a non-directed random 
walk through the accessible configuration space. 

which point only small fluctuations from equilibrium are 
permitted. 

However, no constraint is present for our entangled 
state, which is free to execute a random walk and roam 
over the accessible region of parameter space 72, for ex- 
ample a given initial configuration Q may start in region 
V, but is free to move towards and away from equilib- 
rium and can visit every point in the accessible parameter 
space, see Fig. [5l 

A global pure state 

We now present a highly entangled pure state for which 
the thermodynamic arrow is removed for interactions 
between any two of its subsystems, and so allows the 
most stark deviation from standard thermodynamical be- 
haviour. 

In general for N qubits in locally thermal states there 
exist a vast number of consistent pure states. However, it 
is straightforward to see that any set of thermal marginals 
{pi\iL\ can De obtained from a global pure state with a 
fixed energy E, of the form 

|*) = ^ Xi Xi\0 ■ ■ ■ 0) + y/1 - E/(N - 1)|1 • • • 1) (8) 

i 

with Xi the Pauli operator on qubit z, and each xi real 
and obeying the condition J2i x i — E/(N — 1). 

Furthermore, these parameters are related to the pa- 
rameters for Vn through A^ = Y^j^i x<2 j- 

The parameters X{ specify a point on a hypersphere 
in W N of radius yjEj (N — 1) and any unitary transfor- 
mation that satisfies the energy constaint Qi = 
then corresponds to O(N) transformations on the sub- 
space spanned by {X{\0 • • • 0)}fL 1 . Heat transfer be- 
tween any two qubits a and b can be mediated by an 
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FIG. 6: Ext angled system: The heat flow pattern into 
A,B,C for the correlated state pabc with Xa = 0.15, Xb = 
0.2, Ac = 0.3, 7 = 0.4 and parameters t for the interaction 
Vab and s for Vbc- Here a light-coloured plane is shown 
through Q — to highlight those regions where heat leaves 
the subsystem Q < 0, while the dark regions visible above 
this plane correspond to heat entering the system Q > 0. 



energy-conserving interaction Hamiltonian such as V a b = 
(YaXb — X a Yb)/2 and the evolution is then described by 
the unitary \J a b = exp[— iO/2(Y a X} ) — X a Yb)]. This uni- 
tary allows heat flow in either direction between the ther- 
mal reduced states p a and pi>. 

While the reduced states p a -> pb are thermal, it is not 
the case that p a b = Tr^ a ^|\I/)(\I/| is a direct product of 
thermal states. The 2-qubit process we gave, detects the 
quantum correlations present between a and b that are 
not accounted for by merely probing the qubits individu- 
ally. A state such as (|5J) could clearly be used by a global 
demon to reverse standard thermodynamic behaviour for 
a local demon, as described in the introduction, with an 
interaction such as V a b causing heat to flow both into and 
out of a thermal state. In addition, we see that quantum 
correlations in a pure state allow far more freedom than 
would be allowed by classically correlated states. 



A candidate experimental state? 

We now consider an interesting system consisting of 
3 qubits A, B, and C in a line, where A is allowed to 
interact with B and B is allowed to interact with C. The 
state of the total system is given by a mixed state Pabc-, 
and has the property that pab — Pa(Ta) ® Pb(Tb) and 
Pbc = Pb(T b ) ® pc(T c ) for thermal marginals pi(Ti), 
i = A,B,C with T A <T B < T c . Hence the pairs AB 
and BC obey the thermodynamic arrow individually and 
we would expect heat only to flow from C to B and from 
B to A. 
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FIG. 7: Uncorrelated system: The heat flow pattern into 

A, B, C for the ordinary product state pA® pB pc, with pa- 
rameters Xa — 0.15, Xb = 0.2, Ac = 0.3, t for the interaction 
Hab and s for ii^c- Here a light- coloured plane is shown 
through Q — to highlight those regions where heat leaves 
the subsystem Q < 0, while the dark regions visible above 
this plane correspond to heat entering the system Q > 0. 

However, if pac is given by 

Pac = i((7 + A c -A A )|10)(10| + ( 7 -A c + A A )|01)(01| 

+ V7 2 -(Ac-A A ) 2 (|10)(01| + |01)(10|) 

+(Aa + a c - 7)|00)<00| + (2 - x A - x c - t)|h)(h| 

+ V(Aa + A c - 7)(2 - A A - A c - 7)(|00)(11| + |11)(00|)) 

then we are able to exploit correlations between A and C 
to for example transfer heat from A to B. The state p^c 
at first glance looks complicated, but it may be obtained 
from a 3 qubit state \^}acd of the form (|5J) by tracing 
out D; the parameters A a, Ac, 7 must then obey similar 
relations to (0. 

The total state now takes the form pabc — Pac ® Pb 
and we assume that the interaction Hamiltonians are 
once again given by Vab — \{X aYb — YaXb) and 
Vbc = \{XbYc — YbXc). In the event that only one of 
these interactions is switched on we have that the ther- 
modynamic arrow is obeyed and heat always flows in the 
direction C — )• B — )• A. However when both interactions 
are on, the heat flow becomes more complicated. 

If we assume the system pabc evolves under the uni- 
tary U(t,s) = exp[— UVab — isVBc], then Fig. [6] and 
Fig.[?lshow the contrast between heat flows for the entan- 
gled state pabc and for the uncorrelated product state 

Pa® Pb ® pc- 

For s = heat can only flow from B to A, however as 
can be seen from Fig. [6] by switching on s we are able 
to have system A, the coldest subsystem, emit heat into 

B. Fig. [71 shows the expected heat flow patterns given 
the standard uncorrelated product thermal states. Here 




FIG. 8: The first graph shows the difference in heat pattern 
for A in the two cases of an entangled state and a separable 
state. Here Q P A is the heat gained by A in the product state 
case, while Q e A is the heat gained by A in the entangled state 
case. Once again, light regions correspond to negative values, 
while dark regions correspond to positive values. The second 
plot shows the regions of s-t parameter space in which A has 
lost heat and both B and C have gained heat. 

the heat flow into A is always positive (dark-coloured 
regions), while the heat flow into C is always negative 
(light-coloured regions). It is clear that the presence of 
entanglement can dramatically alter this pattern. Fig. [8] 
shows the difference AQa = Q P a ~ Qa between the flow 
for A in the two cases of a product state (Q^) and an 
entangled state Q e A . It is evident from this that the en- 
tanglement can increase the heat into A as well as out 
of A. Fig. [8] also shows the region of parameter space 
for which heat leaves A and both B and C gain heat, 
which corresponds to a clear violation of the expected 
thermodynamic behaviour. 

We can also ask if the system can act as an entangle- 
ment witness under heat flow. By analysing the mutual 
information for the system we can show that range exists 
in the parameters (A^,Ac,7) for which the correlations 
present in the state pabc exceed log 2 and so allows heat 
flows that cannot be attributed to classical correlations 
between A and C. Fig. [9] shows this region as part of the 
diamond polygon parameter space discussed earlier. The 
set of states for which it may act as an entanglement wit- 
ness is shown as the semi-transparent light-coloured re- 
gion. However, it must be emphasized that for a generic, 
correlated mixed state it is not always possible to uni- 
tarily eliminate all of the correlations present as the set 
of unitary orbits containing product states is a low di- 
mensional sub-manifold of the manifold of all unitary 
orbits [3 lj. 

A quiver of thermodynamical arrows 

As already mentioned in general there exists a multi- 
tude of global states p tot associated with a given combi- 
nation of thermal marginal states. For a given collection 
of marginals the full state can have a vast array of dif- 
ferent consistent correlations. The set of such consistent 
states Q contains states that differ in the correlations 
that exist between different combinations of subsystems. 
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FIG. 9: The region of parameter space for which reversal of 
the thermodynamic arrow in pabc can act as an entanglement 
witness is shown in semi-transparent light region, while the 
remainder of the space only possesses correlations that could 
be attributed to a classically correlated state. 

Consequently, given subsystems Si , . . . S& , the reduced 
state p\...k can vary for different global states p tot in Q. 

It is easy to see that while we might demand the 
individual subsystems Si, ... , Sn to be thermal states, 
this thermality may appear at various levels of coarse 
graining on the collection of subsystems. For example 
we might have for a pair of subsystems S a and S5 that 
p a b = Tr[p tot ] = p a (8> phi however for a triple of sub- 
systems S a , S&, S c the state p a b c is entangled. For such a 
situation, the thermodynamic arrow is present for any en- 
ergy transfer process involving only the two subsystems a 
and 6, however for processes involving three or more sub- 
systems the correlations present allow violations of the 
thermodynamic arrow. 

Consequently, the different possible correlation struc- 
tures that can occur in a state p tot G Q correspond to 
a hierarchy of thermodynamic arrows. The set of states 
Q can be partitioned up as Q = Ai U A2 U • • • , where 
a state p is in Ak if there exists subsystems S ai , . . . , S ak 
such that 

k 

E^>° (9) 

is observed for all transformations on p ai ...a k and no 
larger set of subsystems exist with this property. It is 
not clear, however, if a non-trivial partitioning can be 
defined in which, for a state in a given class, the thermo- 
dynamic arrow holds for all marginals on k subsystems. 

This classification of states in terms of their largest 
thermal marginal depends on the particular correlation 
structure of the global state without reference to where 
the subsystem are situated, however in more practical 
situations it makes sense to include a notion of locality. 
For example, we might consider a collection of subsys- 
tems Si,S2, . . • located at various points in space. For 
a given subsystem S a we may define the notion of an 
"arrow range" , R a , which is the largest r a such that the 



reduced state p(r a ) on all subsystems within a distance 
r a of S a is a thermal marginal. For the entire state p tot we 
may then simply define a characteristic arrow range R as 
the average over the set {Ri} of ranges for each subsystem 
Si. States with large R do possess correlations, however 
these correlations are difficult to access in practice, re- 
quiring many coordinated, local pairwise interactions Vij 
in order to generate an effective multipartite interaction 
over the correlated state. On the other hand, states with 
small values for R, for example states of the form (|SJ) or 
the state pabc considered in the previous section, have 
correlations that are more easily accessible through local 
pairwise interactions. This provides a simple generaliza- 
tion of the 'local' Maxwell demon introduced in the first 
section. A demon should not only be finite in terms of 
its memory resources, but also in terms of the correlation 
range R that it can probe. 

CONCLUSION 

Time, with its inexorable flow, is one of our oldest mys- 
teries. In stark contrast, it is only in the last century that 
we have become aware of the quantum mechanical prop- 
erties of Nature, and only in the past few decades that 
some of the subtleties and power of quantum correlations 
have revealed themselves. 

In this article we have examined some aspects of entan- 
glement correlations, how they can reverse the standard 
thermodynamic arrow of time, and have shown that in 
this regard thermodynamic transformations may act as 
entanglement witnesses. 

We analysed how highly entangled multipartite states 
can make ordinarily forbidden thermodynamic states ac- 
cessible. Different correlation structures can give rise to 
the same set of local thermal states, and so processes 
involving several subsystems are in general required to 
exploit the entanglement present. The set of transforma- 
tions on a system then possess a hierarchy of thermody- 
namic arrows, whether considered globally or in terms of 
local correlation ranges. The key aspect is to allow inter- 
actions that activate the correlations present, which we 
demonstrated with a mixed state example pabc in which 
the thermodynamic arrrow is in place for interactions in- 
volving only A and B or for B and C, but by switching 
on an interaction with system C this can be dramatically 
modified. 

It would be of interest to further explore the multi- 
partite correlation structures that can occur in a quan- 
tum state - for example we might speculate that in the 
early, dense universe that R is small enough to allow any 
random physical interactions to exploit the correlations 
present, producing a gradual disappearance of the ther- 
modynamic arrow the closer we get to the initial state of 
the universe. Of course such notions of locality rely on 
a classical spacetime background, and so we could only 
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push the issue of the special initial conditions into the 
sub-Planck scales, where a complete theory of quantum 
gravity would be required. 

Another setting of interest is that of quantum field the- 
ory, where thermality arises naturally due to observers 
possessing causal horizons, such as in the case of black 
holes or for accelerated observers. Once again, thermal- 
ity can be viewed as arising from entanglement across 
the horizon, however it is unclear if these correlations 
could actually be utilized in any thermodynamic trans- 
formations. At a deeper level, it may also be fruitful 
to explore connections with the thermal time hypothesis 
due to Connes and Rovelli [32[ or recent work by Pad- 



manabhan [33| that attempts to relate the cosmological 
arrow of time for our expanding universe with the ther- 
modynamic arrow by exploiting the thermodynamics of 
local spacetime horizons. 

Finally, in a more grounded setting, it would be of 
interest to see if a mixed state Pabc, ot one like it, could 
be realised experimentally to see how easily one might 
exploit quantum correlations present in a state to affect 
its thermodynamic behaviour. 
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